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Electrical current fluctuations in a single-channel quantum point contact can produce photons 
(at frequency uj close to the applied voltage V x e/h) which inherit the sub-Poissonian statistics 
of the electrons. We extend the existing zero-temperature theory of the photostatistics to nonzero 
temperature T. The Fano factor J" (the ratio of the variance and the average photocount) is < 1 for 
T < Tc (antibunched photons) and > 1 for T > Tc (bunched photons). The crossover temperature 
Tc ~ Acj X h/ks is set by the band width Aw of the detector, even if hAuj <C eV. This implies 
that narrow-band detection of photon antibunching is hindered by thermal fluctuations even in the 
low-temperature regime where thermal electron noise is negligible relative to shot noise. 

PACS numbers: 73.50.Td, 42.50.Ar, 42.50.Lc, 73.23.-b 



I. INTRODUCTION 

In classical mechanics, time dependent fluctuations of 
an electrical current produce photons with the Poisson 
statistics of classical particles.'^ The variance Varn of 
the number n of photons detected in a time t^ct is then 
equal to the mean (n). Quantum mechanics changes 
the photostatistics.'^] The bosonic nature of the photons 
would naturally lead to photon bunching, with Var n > 
(n). Photon antibunching, with Var < (n), is also possi- 
ble, if the photons can inherit the sub-Poissonian statis- 
tics of the electrons.!^ 

It is an experimental challenge to detect antibunched 
photons produced by electronic shot noise in a quantum 
conductor.^ The theoretical prediction^ is that pho- 
tons emitted by a single-channel quantum point contact 
should have a Fano factor = Vain/{n) smaller than 
unity at zero temperature, for frequencies uj close to the 
applied voltage V x e/h. More specifically. 



T= 1 



(7oAc^)r(l - r) 



(1.1) 



for photodetection with efficiency 70 in the frequency in- 
terval [eV/h—Auj, eV/h). The transmission probability r 
through the quantum point contact is assumed to be en- 
ergy independent on the scale of eV. Eq. ( 1.1 1 is derived 
in the limit of weak coupling (7oAaj ^ 1) of electrons 
to photons, so that the deviations from Poisson statis- 
tics remain small. It is also assumed that the photons 
can be detected individually, see Ref. [7|for an alternative 
detection scheme. 

It is the purpose of the present paper to extend the the- 
ory of Ref.|n|to nonzero temperatures, in order to identify 
the conditions on the temperature needed to observe the 
photon antibunching. Clearly, photon bunching should 
take over when the electrical shot noise drops below the 
thermal noise, which happens when fc^T becomes larger 
than eV. While ksT < eV is the condition for photon 
antibunching in the case of wide-band detection, a more 
stringent condition fc^T < hAuj holds for narrow-band 
detection. 



More precisely, we obtain a crossover temperature 
Tc w hAuj/AkB at which J" = 1 for Aw <C eV . In this 
low-temperature regime shot noise still dominates over 
thermal noise, yet the photon antibunching is lost. One 
qualitative way to understand this is, is to compare the 
coherence time tcoh — l/Acj of the detected radiation 
with the coherence time It — h/kBT of thermally ex- 
cited electron-hole pairs. For tcoh > tr the detected pho- 
tons result from many uncorrelated electron-hole recom- 
bination events, and the one-to-one relationship between 
electron and photon statistics is lost. 

In the next section we give the nonzero temperature 
generalization of the theory of Ref. and then in Sec. 
|III| we specialize to the shot-noise regime fc^T ^ eV . 
General results in both the shot noise and thermal noise 
regimes are presented in Sec. |IV[ Technical details are 
summarized in the Appendices. 



II. GENERATING FUNCTION AT NONZERO 
TEMPERATURE 



We seek the nonzero-temperature generalization of the 
formulcP 



N 



m— 1 



J] Det l + T„Je^e^ -1 



(2.1) 



for the factorial- moment generating function F{£) of the 
photocount. We first introduce the notation and then 
present the required generalization. 

The photons are produced by time-dependent current 
fluctuations in a quantum point contact, characterized by 
transmission eigenvalues Ti, T2, . . . Tat, with N the num- 
ber of propagating electronic modes. The current flows 
between two reservoirs, with Fermi functions 

fL{e) = (1 + exp[(£ -eV- EF)/kBT]y\ (2.2) 
fR{e) = (1 + exp[(e - EF)/kBT]Y\ (2.3) 
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The current fluctuations can be due to thermal noise (at 
temperature T) or due to shot noise (at a voltage V 
applied over the point contact). We take the transmis- 
sion eigenvalues Tm as energy independent in the range 
max(ey, ksT) near the Fermi energy Ep. 

The photons are detected during a time idot in a narrow 
frequency interval Auj around frequency 57, as determined 
by the detection efficiency 7(aj). Antibunching of the 
photons requires that ft is tuned to the applied voltage, 
n ~ eV/h. (In the following we set h and e equal to 
unity.) 

The average (• • • ) in Eq. (2.1 1 indicates a Gaussian 



integration over the complex numbers Zp, 



n 



7p 



(Pzr,e 



-7p|2pl 



(2.4) 



The matrix Z has elements 



pp 



^1/2, 



'■p—p'^p—p' J 



de- 



pending only on the difference of the indices p and 
p' . This difference represents the discretized frequency 



ijp-p' = {p—p')x 27r/tdet of a photon emitted by an elec- 
tronic transition from energy Sp to Spi and detected with 
efficiency jp-p' — (27r/t(jet)7('^p-p')- Since j{uj) = for 
w < 0, the matrix Z is a lower-triangular matrix. The 
discretization of frequency and energy is eliminated at 
the end of the calculation, by taking the limit tdct ^ oo. 
The expansion 



GO 

k=0 



(2.5) 



of F{^) in powers of ^ gives the factorial moments 
(n'^) / — {n{n — l){n ~ 2) ■ ■ ■ {n — k + 1)) of the num- 
ber of detected photons. Antibunching means that the 
variance of the photocount Var n — {n?) — (n)^ is smaller 
than the average, or equivalently that the Fano factor 
T ^ Varri/(n) < 1. 

As outlined in App.jAj at nonzero temperature we have 
instead of Eq. (2.1) the generating function 



N 



Y[ Det 



\m—l 




(2.6) 



The Fermi function fhi^) in the left electronic reservoir 
is contained in the diagonal matrix /i, with elements 
Ul)pp' = Spp-hisp), Ep = p X 27r/tdot. Similarly, the 
Fermi function fR{s) in the right reservoir is contained 
in the diagonal matrix //j. 

Following the steps in App.|Aj the expression (2.6 1 can 
be reduced to the more compact form 



N 



(2.7) 

R.M = 

The zero-temperature limit (2.1 1 follows from 



with the definitions = ^ ^ II, Ir = ^ ^ 



E 



Eq. (2.7 1 by setting /l = 1, = in the energy interval 
< e < Ep + V . (There are no current fluctuations 



outside of this energy interval for T — 0.) 



III. SHOT NOISE REGIME 



The result (2.7) holds for any temperature, provided 



that the energy dependence of the transmission eigenval- 
ues may be neglected. In particular, it describes both 
thermal noise and shot noise. A simpler formula is ob- 
tained in the shot noise regime ksT <^ V. Thermal 
noise can then be neglected and only the finite temper- 
ature effects on the shot noise are retained. We assume 
Aw <C — V^, so even if ksT ^ V, the relative magni- 



tude of Auj and ksT is still arbitrary. 



A. Generating function 



The first simplification in this regime is that we may 
set fRe-^fL 0, since /fl,(£)/L(e') ^ for e' < e. Eq. 
(2.7) reduces to 



N 



FiO - ( n Det e-^' + T„jLMfi 



(3.1) 



where we have multiplied by Det e^^^ — 1. 

The second simplification is that we can ignore energies 
separated by pV with p > 2, because V is the largest 
energy scale in the problem. Since and Z^''' connect 
energies separated by pQ ~ pV, we may set ZP,Z''^ — >■ 
for p>2. From Eq. (3.1) we arrive at 



N 



FiO = { n Det l-Zt+r„Ji(Z + Zt)/^ 



(3.2) 



Following the steps in App. [B] the determinant may be 
rewritten in the more convenient form (bilinear in Z, Z^), 



N 



FiO = ( n Det 1 + T^(l - T^)fLZfRZ^ 



(3.3) 
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B. Moment expansion 



C. Crossover from antibunching to bunching 



The generating function (3.3) is of the form F{^) ~ 
Det (1 + Xfn) with of order An expansion in 
powers of ^ can be obtained by starting from the identity 

n^Det (1 + X„0 = exp[E„Tr ln(l + X„,)] , (3.4) 

and expanding in turn, the logarithm and the exponen- 
tial Up to second order in ^ we have the expansion 



F{0 = 1 + (E™TrX,„) - i(E™TrX2,) 
+ K(E™TrX™)')+0(e3), 



(3.5) 



from which we can extract the first two factorial mo- 
ments, 

F{0 - 1 + ^{n) + '^eiin') ~ (n)) + ©(C^). (3.6) 

We perform the Gaussian averages and obtain the aver- 
age photocount (n) and the variance Varn — (n^) — (n)^ 
in the shot noise regime, 



/ \ ^dct 

(n) = 

^ ' 27r 

Var n = (n) 



Si I du^iiu) I defL{e + Lo)fR{e), (3.7) 

2 



2tt ^ ' 



tdct 
^dct 



de 
de 



7(c^) j defL{e + Lo)jR{e) 

■ 2 

/L(e) / dtJ7(w)/ij(£-a;) 
/fl(e) / da;7(a;)/L(£ + w) 



We have defined 



(3.8) 



(3.9) 



Since the two reservoirs are at the same temperature, 
we can write /l(£) = f{e — V — Ep) and fn = /{Ep — s) 
in terms of a single Fermi function 



/(e) = (l + e^/'=-^)-i. 



(3.10) 



We abbreviate r {e,u j) — j{uj)f{e)f{uj — e — V) and can 
then write Eqs. (3.7) and (|3.8|) in the compact form 



^dct 

2^ 



Si J duj J deT{e,u}) 



(3.11) 



Var n 



tdct 

27 



duj / deT{s,u!) 



S^ / de'r{e',uj)~2S2 / duj'T{e,uj') 



(3.12) 

The difference Var n — (n) contains a positive term cx Sf 
and a negative term cx 5*2. The sign of this difference 
determines whether there is bunching or antibunching of 
the detected photons. 



To investigate the crossover from antibunching to 
bunching with increasing temperature, we take a block- 
shaped response function 



7(0;) 



70 if V - Aw < u! <V, 
otherwise. 



(3.13) 



In the low-temperature regime ksT Aw the function 
r(e, w) then has a block shape as well and we recover the 
results 



{n} = ^o^^2 ' 



Var n — {n) 



'^{loA.rl{Sl~2S,) 



of Ref. |61 These correspond to a Fano factor 
^ = 1 + |7oAw(5i - 2S2/S1). 



(3.14) 
(3.15) 

(3.16) 



For a single-channel conductor S2 = Sf, so there is anti- 
bunching (J^ < 1) at low temperatures. 

At high temperatures ksT ^ Aw, but still in the shot- 
noise regime k^T <^ V, we may substitute T{e,Lu) 
~-f{uj)kBTdf{e)/de into Eqs. ( |3.11| ) and ( |3.12| , which 
gives 



tdct Aw 
zvr 

Var n — {n) = 
The Fano factor 



jnksTSi, 



tdct Aw 
27r 



(jokBTySf. 



F = 1+ joksTSi 



(3.17) 
(3.18) 

(3.19) 



is now > 1 — hence there is photon bunching. 

The crossover temperature Tc, at which J- = 1, can be 



calculated numerically from Eqs. (3.11) and (3.12). In 
the single-channel case, when 52 = Sf, we find 



ksTc « 0.25 Aw. 
The crossover is shown graphically in Fig. [l] 



(3.20) 



IV. BEYOND THE SHOT NOISE REGIME 

In the previous section we assumed ksT <^ V (shot 
noise regime). For arbitrary relative magnitude of ksT 
and V, the general formula (2.7) can be used. With 



the help of Eq. (3.4), this general expression of the form 



Det (1 -I- X) was expanded to second order in powers of 
In this case however, since X ~ 0{y/^), terms up to 
order X'^ had to be retained. The first two moments of n 
are obtained as integrals over energy and frequency, sim- 
ilar to Eqs. (3.111 and (3.12) but containing many more 
terms in the integrands. The results shown in Figs. [2] and 
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FIG. 1; Crossover with increasing temperature from anti- 
bunching (Fano factor < 1) to bunching (T > 1) of the 
photons produced by a single-channel quantum point contact 
in the shot noise regime (fesT <^ V). The solid curve is calcu- 
lated from Eqs. (3.11 1-( 3.13 1. The dashed line is the asymp- 



tote (3.191. The crossover temperature Tc from Eq. (3.201 is 



indicated. 




kBT//\uj 



FIG. 2: Same as Fig. [1] but now without making the restric- 
tion to the shot noise regime (so without assu ming fesT <C V^). 
The two solid curves are calculated from Eq. (2.71 for two val- 



ues of Aoj/fcaT (both for the single-channel case with trans- 
mission probability r = 0.5). Both curves converge to the 
shot noise result at low temperatures (shown dashed). 



[3] are for the case iV = 1, Ti = r of a single channel, and 
for the box-shaped response function (3.131. 



As expected, all curves converge to the shot-noise re- 
sults when kgT ^ V (shown dashed). At higher tem- 
peratures, the Fano factor lies above the shot noise limit 
due to the appearance of thermal noise. The tempera- 
ture Tc at which antibunching crosses over into bunch- 
ing, so when = 1, follows the shot noise limit (3.20) for 
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FIG. 3: Dependence of the crossover temperature Tc (at which 
= 1) on the band width Au). The points are calculated from 
Eq. (2.7 1 for three values of the single-channel transmission 
probability r. For Ao; <C V all points converge to the shot 



noise limit (3.201 (dashed line) 



narrow-band detection (Aw <C V) . With increasing band 



width, Tc drops below the shot noise limit, in particular 
for small transmission probability r. For r = 0.5 the shot 
noise limit remains accurate even for band widths Acj as 
large as V/2. 



V. CONCLUSION 

In conclusion, we have investigated the effects of a 
nonzero temperature on the degree of antibunching of 
photons produced by current fluctuations in a quantum 
point contact. Antibunching crosses over into bunching 
as a result of thermal noise in the point contact, but 
this is not the dominant effect in the case of narrow- 
band detection. In that case the finite coherence time 
of electron-hole pairs governs the transition from pho- 
ton antibunching to photon bunching, which occurs at a 
temperature ksTc — Aw even if ksTc V (so even if 
thermal noise is negligible relative to shot noise). 

The optimal conditions for the observation of anti- 
bunched photons are reached for a band width Aw « V/2 
and a transmission probability r « 1/2 through a single- 
channel quantum point contact. In that case /csTc « V/S 
has the largest value at any given applied voltage. 
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Appendix A: Derivation of the generating function 
at nonzero temperature 

We briefly describe how the analysis of Ref. [S] can be 
generalized to nonzero temperatures, in order to arrive 
at Eq. (2.6). Referring to the equations in that paper, 



the first equation which changes is Eq. (5), which now 
reads 



m = { 



a^DZa b'^DZb b^DZ''b -a'' DZ'' a\ 

e e e e ). 



(Al) 



We have Bij = Sijfi, with /, the Fermi occupation 
number in channel i. The matrix A is given by A = 
e-^e^e^^e^' - 1, with X ^ -DZ and Y = S^DZS. 
Notice that XP = D{-Z)p and YP = S^DZPS. 

We now make the assumption of an energy independent 
scattering matrix, so S, commute with Z, Z^ . The 
determinant is invariant under a change of basis, and by 
working in the eigenchannel basis we can reduce S* to a 
2x2 matrix Sm for each eigenchannel. 



The four factors correspond to the four current operators 



that need to be taken into account: I- 



in' ^outJ ^out 



/out, and 



The operator creates an incoming electron, while 6^ 
creates an outgoing electron. The matrix D projects on 
the right lead, where the current is evaluated. (Since D 
commutes with Z, we can write DZ instead of DZD.) 
One can relate b = Sa, with S the unitary scattering 
matrix, so one can write the entire generating function 
in terms of the operators a. The expectation value (• ■ ■ ) 
is both an expectation value over the fermion operators a, 
as well as the average over the Gaussian variables Z, Z^ . 

Following the steps of Ref. i6j we calculate the expec- 
tation value of the fermion operators by means of the 
identity 



JJe"^-^"") =Det(l + AB), 



(n. 



1, Bij = {a\ai 



(A2) 
(A3) 



V J- ?T). 



(A4) 



1,2,...N the transmission eigenvalue. 



with Tm, m 

The matrix structure of /, D, and Z in this basis is 



/ = 



h 
Jr 



D = 




1 



z = 



Z 
z 



(A5) 



Substitution of Eqs. (|A2|-(|A5| into Eq. (AA\ le ads after 
some algebraic manipulations to the result (2.6 1. 



The determinant in Eq. ( 2.6 1 can be reduced by means 
of the folding identity 



Det 



Mil Mi; 



^M2i M22 
leading to 



DetAfii Det (M22 - M2iAff/Mi2), 

(A6) 



N 



F{0 = { n Dct[l +T,„/i(e^e^' - l)]Dct( 1 + r„/«(e-^e-^' - 1) 



We continue the reduction of the determinant, using first the identity 



1-1 ^ /'„-Z+ „Z\ 



e — e 



then multiplying the determinant by Dete'^* = 1, and finally combining the product of three determinants into a 
single determinant. In this way we eliminate the matrix inversion, arriving at 

m = ( n [1 + TMe-^e-^' 1)] e^' [l + T„J, [e^ e^' 1)] 

-Z „ZUf („ZZ'^ 



(A7) 



(A8) 



+ Tmil - T„,)Me-^ - )fL{e^e^ - 1) 



N 



[] Det 1 + Tm[{l - fR)e^ h ~ fne-^l - /l)]( 



(A9) 



This is Eq. (2.7) in the main text. 



Appendix B: Derivation of the generating function the steps required to arrive at the bilinear form (3.3 I 
in the shot noise regime 



Starting from the expression (3.2) for the generating 
function in the shot noise regime /c^T <C we give 
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We group terms with Z and with Z"^ in the matrices 
Am = TmfLZfR and = T^fLZ^h - Z\ so that Eq. 



(3.2) can be written as 

AT 



(Bl) 



Eq. (3.3) fohows by noting that ZfRf^ — Zfn for 
fcsT ^ 17 ~ 1^, since the Fermi function in this term 
is evaluated at energies near Ep , where it can be replaced 
by unity. Similarly Z'^ JnfL Z'^ Jl-, since Jr is evalu- 
ated at energies near Ep + V where it can be replaced 
by unity. 



Because energies separated by with p > 2 can be 
discarded, we may set — 0, — 0. For any pair of 
matrices A, B which square to zero, one has the identity 



Det {1 + A + B) = Dct(l - AB). 
This leads to 

/ ^ 

FiO = { n Det(l + T„,ZJrZUl 



fRfhZ^ fRfh] 



(B2) 



(B3) 
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